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Purcell's scallop theorem defines the type of motions of a solid body - reciprocal motions - which 
cannot propel the body in a viscous fluid with zero Reynolds number. For example, the napping of 
a wing is reciprocal and, as was recently shown, can lead to directed motion only if its frequency 
Reynolds number, Re/, is above a critical value of order one. Using elementary examples, we show 
the existence of oscillatory reciprocal motions which are effective for all arbitrarily small values of 
the frequency Reynolds number and induce net velocities scaling as Re" (a > 0). This demonstrates 
a continuous breakdown of the scallop theorem with inertia. 

A large variety of biological movements occur in a fluid environment, from swimming bacteria to whales. In many 
cases, the study of fluid forces is crucial to the understanding of animal locomotion [l[ El, E[ 0, El Q ■ Because of the 
large range of relevant length scales in biological motility - eight orders of magnitude in size, from less than a hundred 
nanometers to tens of meters - fluid mechanics occurs in distinct regimes with important mechanical consequences. 
On small length scales, the relevant Reynolds number is usually very small (Re 10~ 4 for swimming E. coll) and 
viscous forces are dominant. This is the Stokesian realm of swimming microorganisms such as bacteria, spermatozoa, 
and ciliated cells. At the opposite end of the range of length scales, the Reynolds numbers are typically very large 
(Re 10 7 for a swimming tuna) and inertial forces are dominant. This is the Eulerian realm of flying birds and 
swimming fishes. In this paper, we address the transition from the Stokesian to the Eulerian realm, and show that, 
in some situations, this transition can take place continuously with an increase of the relevant Reynolds number. 

In his 1977 lecture Life at low Reynolds numbers, Edward Purcell introduced the "scallop theorem" Q. He observed 
that the Stokes equations, which govern fluid flows at zero Reynolds numbers and are both linear and independent 
of time, are identical under time reversal. Consequently, there exists a certain geometrical class of motion (or, more 
generally, actuation of a solid body), termed "reciprocal motion", which cannot lead to any locomotion in this limit. 
A reciprocal motion (or actuation) is a motion in which the geometrical paths followed by various material points on 
the body are identical when viewed under time reversal. By symmetry, such motion can only lead to a net movement 
equal to minus itself, and therefore, no net movement at all (see also Refs. The simplest example of a reciprocal 

motion is a periodic motion composed of two distinct parts. In the first part, the body moves in a certain prescribed 
way, and in the second part, the body moves in a manner which is identical to the first-part as seen under time 
reversal. A scallop opening and closing belongs to this subclass of reciprocal motion and, independently of the rate 
of opening and closing, the scallop cannot move. 

Another example of reciprocal motion - or, in this case, reciprocal actuation - is a flapping body. Consider a solid 
body oscillated up and down in translation in a prescribed manner by an external means. Since the motion going up 
is the time-reversal symmetry of the motion going down, the flapping body does not move on average in the limit 
of zero Reynolds numbers. However, large animals such as birds use flapping wings for locomotion, and so clearly a 
thin flapping body must be effective in the Eulerian realm. The question then arises: When does a flapping body - or 
more generally a reciprocal motion - become effective? How much inertial force is necessary to break the constraints 
of the scallop theorem? 

This question was first formulated and studied by Childress and Dudley Q. The mollusc Clione antarctica was 
observed to possess two modes of locomotion. The first is non-reciprocal and uses cilia distributed along the body of 
the mollusc. The second is reciprocal and consists of two flapping wings. The flapping-wing mode was observed to be 



* Email: lauga@mit.edu 



2 



predominant for the large swimming velocities. Using both experimental observations and fluid mechanics models, the 
authors postulated that reciprocal motions are ineffective in producing any net motion unless the relevant frequency 
- or "flapping" - Reynolds number, Re/, is sufficiently large (order unity). In other words, the transition from no- 
motion to motion occurs at a finite value of Re/ and the breakdown of the scallop theorem is discontinuous. This idea 
was subsequently studied in laboratory experiments 0, 0] and numerical simulations ~U^, 11 1 of flapping symmetric 



bodies, both of which confirmed the transition to directed motion as a symmetry-breaking instability occurring at a 
finite value of the frequency Reynolds number, as well as the robustness of this transition to a change in a variety of 
geometrical and mechanical parameters. 

In this paper, we consider a series of elementary oscillatory reciprocal motions of a solid body with broken spatial 
symmetries and show that they become effective in producing a net translation of the body for arbitrarily small 
values of the frequency Reynolds number, with induced velocities scaling as Re/ (a > 0, inertia! creep [l^). This 
demonstrates a continuous breakdown of the scallop theorem with inertia. 

The examples we propose rely on classical results of lift forces for the motion of spherical particles at small Reynolds 
number 0,0,0, 1 1 51 ] - We consider a solid spherical particle (density p pi radius a) oscillating with frequency u> and 
amplitude d in a fluid of density p and shear viscosity p. The three different setups we propose are described below, 
and we start by some general remarks. In the case of purely translational motion, including the effect of inertia on the 
particle motion can be done in a number of limits, as there arc in general three relevant Reynolds numbers. Firstly, 
the unsteady term in the Navier-Stokes equations scales as puiUo (where Uq = doj is the typical speed of translation), 
and is smaller than the typical viscous term, of order pUo/a 2 , by a factor of Re w = a 2 uj/v, where v = p/p is the 
kinematic viscosity. Secondly, the nonlinear advective term in the Navier-Stokes equations scales with pU^/a, and is 
smaller than the viscous term by a factor of Re/ = aUa/v — aduj/v, which is the flapping (or "frequency") Reynolds 
number @, Thirdly, the particle inertia is quantified by a particle Reynolds number, Re p = p p a 2 uj/p, the ratio 

of the typical rate of change of the particle momentum, p p a 3 Uou>, to the typical viscous forces on the particle, pallo. 
In this paper, we will consider the asymptotic limit where 

{Re p ,Re w } < Re/ < 1, (1) 

so that the motion of the flapper is quasi-static and the leading-order departure of the fluid forces from the Stokes laws 
is due to the nonlinear advective term in the Navier-Stokes equations [28|. The limit described by Eq. ([1]) is equivalent 
to that of small frequency Reynolds number (Re/ <C 1) and larg e flapping amplitude (a/d <C 1 and a/d <C p/p p )- 



and larg 

Note that this is a different limit from the work in Refs. [1, [s, njlll| where body inertia likely played an important 



role. We consider below three examples of such large-amplitude low-Re/ reciprocal flapping which lead to directed 
motion for arbitrarily small values of Re / . 

The first example is that of a flapper near a wall. Specifically, we consider the reciprocal oscillation in vertical 
position of the solid sphere with velocity U(t) = U(t)e x parallel to a stationary solid surface and free to move in the 
y and z directions (see notations in Fig.QJi). In the Stokes flow limit (Re/ = 0), the sphere experiences no lift force, 
and remains at a constant distance, h, to the solid surface. The first effect of inertia on this problem, in the limit set 
by Eq. (fll), is the appearance of a lift force, directed away from the solid surface, and independent of the sign of U(t) 
13 , 17. lid]. Such a limit is captured when the Osccn length scale v/Uq, the distance away from the sphere where 
inertial forces become important, is much larger than all relevant length scales of the problem, i.e. the sphere radius, 
a, and its distance to the surface, h. In the simple case where a <C h <C v/Uo, the lift force leads to a low-Reynolds 
number lift velocity for the particle 
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Vx(t) = V ± (t)e y , V± (t) = ^^^, (2) 

always directed away from the surface. For an oscillatory motion, U(t) = [/ocoswi, the lift velocity away from the 
surface averages over one period to 

M = -Re/. (3) 
U D 64 ' K ' 
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FIG. 1: Three examples of reciprocal forcing leading to translation of a solid body (sphere of radius a) for arbitrarily small 
values of the frequency Reynolds number. In each example, the reciprocal motion is composed of the periodic repetition of two 
distinct parts, with the second part (bottom) being identical to the first part (top) as seen under time reversal: (a) Oscillation 
in vertical position of a sphere parallel to a solid surface leads to motion perpendicular to the surface; (b) In-phase oscillations 
in translation and rotation of a sphere leads to motion perpendicular to both the directions of translation and rotation; (c) 
Oscillation in vertical position of a sphere in an oscillating shear flow (in phase) leads to motion perpendicular to the direction 
of translation; In all cases, V± denotes the (small) sphere velocity induced by inertial forces. In case (a), the distance to the 
solid surface is denoted h. 



A flapper near a wall performing a reciprocal translational motion is therefore able to move forward (away from the 
wall) for arbitrarily small values of the frequency Reynolds number. This inertial migration decreases to zero with 
the first power of the Reynolds number (a = 1), and the Stokes limit is recovered when we formally set Re/ = in 
Eq. ©. 

Our second example is that of a rotating flapper. We consider the case where the solid sphere is oscillating both 
in translation and rotation, with velocity and rotation rates given by U(£) = U(t)e x and S7(t) = fi(i)e z , and is free 
to move in the y and z direction (see Fig. [TJd) . If the two oscillations are in phase, the actuation of the sphere is 
reciprocal, which we will assume here, and no average motion is obtained in the Stokes limit. If £lo is the typical 
magnitude of S~2(£), the rotation Reynolds number Ren = a 2 Q,^/v measures the importance of inertial forces due to 
the rotational motion. In the asymptotic limit set by Eq. ([1]), and for Reo ~ Re/, the first effect of inertia is the 
appearance of a lift force perpendicular to both the directions of translation and rotation 1 
Fl = Tra 3 pfl x U. This results in a low- Reynolds number lift velocity 
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and given by 



/ n xt / \ a 2 U(t)n(t) 
V x (f) = V x (t)e y , V x (t) = y -^±- 
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When U(t) = aSl(t) = Uocoscut, we obtain an average translational velocity, along the y direction, given by 

{V ± ) _ Re/ 
U 12 • 



(5) 



Here again, the reciprocal translational and rotational motion of the solid sphere leads to a directed motion for 
arbitrarily small values of the Reynolds number. The magnitude of this directed motion also decreases to zero with 
the first power of Re/ (a = 1). 

As a final example, we show that these results are also valid when the fluid in the far-field is not quiescent by 
considering a flapper in a shear flow. Specifically, as shown in Fig. [TJ;, we consider the case when the solid sphere 
is oscillating in vertical position with a prescribed velocity, U(t) = U(t)e x , in a shear flow described by the far-field 
undisturbed flow field Uoo = —■y(t)ye x (the center of the sphere is located at y = 0) and is free to move in the y 
and z directions. If the two oscillations are in phase, the motion of the sphere is reciprocal, which we assume here, 
and no average motion is obtained in the limit of zero Reynolds number. We also assume that the sphere is far away 



from the surfaces responsible for the creation of the shear flow and therefore ignore wall effects 13, l2l|, |22| . If 70 
denotes the typical magnitude of j(t), an additional Reynolds number, Re-y = a 2A fo/v, needs to be introduced. Here, 

the first effect of inertia is the a ppe arance of a lift force directed across the undisturbed streamlines [13j, 0,13. The 

1/2 

original study, due to Saffman [23, |24j, calculated this lift force in the limit where Re/ Re^ <C 1, and in this 



case the lift force is moving the sphere in the direction opposite to its translational velocity We consider here the 
same asymptotic limit, together with the limit assumed in Eq. ([T]). In this case, and if U(t) ■ ^(t) > 0, the sphere 
experiences a low-Reynolds number lift velocity given by 

V±(t) = V ± (t)ey, V±(t) = Cl \U(t)\ f^MA) ' ; (6) 



where c\ w 0.343 is a numerical coefficient. For an oscillatory motion U(t) = Uocosuit, and with j(t) = U(t)/a to 
satisfy Saffman's asymptotic limit, we get an average velocity, along the y direction, given by 

( V ±) „ 1/2 

— — = c 2 Re/ , (7) 

where C2 = 2cx j^ 2 (cost) 3 / 2 dt/n w 0.191. As in the previous cases, the actuation of the sphere is reciprocal and yet 
it leads to a directed motion for arbitrarily small values of the frequency Reynolds number Re/. Here, however, the 
magnitude of the induced velocity decreases to zero with the square-root of the Reynolds number (a = 1/2). Also, in 
this case, the motion will continue until the point along the y axis where the local velocity from the shear flow cancels 
out the translational velocity of the sphere. 

As a summary, we have presented elementary examples of oscillatory reciprocal forcing of a solid body leading to 
net translational motion of the body for arbitrarily small values of the frequency Reynolds number, Re/. When the 
frequency Reynolds number is formally set to zero, the effect disappears as dictated by the scallop theorem, but it 
remains non-zero for all non-zero values of Re/. The induced average velocities scale as Re^ (a > 0), corresponding 
to the limit of asymptotically large Strouhal number, St = ujd/ (V±) ~ HcJ a . This demonstrates that the breakdown 
of Purcell's scallop theorem with inertia can take place in a continuous way without a finite onset of translational 
motion. 

As our examples show, a directed motion on the order of the flapping velocity will take place when Re/ ~ 1. 
Moreover, the mechanical efficiencies of the examples above - ratio of the useful work to the total work done by the 
flapper - scale as Re/ Q so that order one efficiencies should also be expected for order one Reynolds numbers. From 
a biological perspective, both these observations suggest that reciprocal gaits are very inefficient for small Reynolds 
number and become advantageous only when Re/ ~ 1. Consequently, and even in the absence of a mathematical 
bifurcation, the onset of an appropriately defined "efficient flapping flight" is expected occur at a finite value of Re / 

Furthermore, it is important to note that all of our examples display some spatial broken symmetries which govern 
the direction of the net motion of the solid body: (a) the location of the wall, (b) the direction of the rotation rate, and 
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(c) the direction of the shear flow. This is somewhat different from the "flapping wing" setup studied experimentally 
in Refs. 0, Q and numerically in Refs. n], 111 where both the shape and the actuation of the wing are symmetric 
and where locomotion is a result of a hydrodynamic instability [29 1. 

Finally, we have considered examples lea ding to net translational motion, but similar examples exploiting lift forces 



and torques on asymmetric particles 
and rotation, of the solid body jiT 
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27] could be devised leading to a net rotation, or combined translation 
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